Carbon Nanotubes as Nanoelectromechanical Systems 
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We theoretically study the interplay between electrical and mechanical properties of suspended, 
doubly clamped carbon nanotubes in which charging effects dominate. In this geometry, the ca- 
pacitance between the nanotube and the gate(s) depends on the distance between them. This 
dependence modifies the usual Coulomb models and we show that it needs to be incorporated to 
capture the physics of the problem correctly. We find that the tube position changes in discrete steps 
every time an electron tunnels onto it. Edges of Coulomb diamonds acquire a (small) curvature. 
We also show that bistability in the tube position occurs and that tunneling of an electron onto the 
tube drastically modifies the quantized eigenmodes of the tube. Experimental verification of these 
predictions is possible in suspended tubes of sub-micron length. 

PACS numbers: 73.63.Nm, 73.23.Hk, 62.25.+g, 46.70.Hg 
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I. INTRODUCTION 

Nanoelectromechanical systems (NEMS) convert elec- 
trical current into mechanical motion on a nanoscale and 
vice versa. They can be viewed as the successorsi of 
microelectromechanical-devices (MEMS) which operate 
at a micron scale and which are found in commercial 
applications. Improved performance is expected from 
NEM-devices due to their small sizes, and higher eigen- 
frequencies. M(N)EMS have already been used for high- 
precision measurements of force 2 «, electric charge^, the 
thermal conductance quantum^, and the Casimir force 5 . 
From a fundamental point of view, NEM-physics is an 
unexplored field in which new phenomena are likely to 
be found. Examples include tunneling through mov- 
ing barriers 6 , additional sources of noisei, and shuttling 
mechanism for transport 8-9-10 . 
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devices made with silicon technology. Carbon nanotubes 
provide an interesting alternative because of their supe- 
rior mechanical properties. They have already been im- 
plemented as nanotweezersiii^, as switches in a random 
access memory devicei^, or as nanoscale actuators^. In 
addition, recent theoretical calculations show that car- 
bon nanotubes can also be used as nanoelectromechani- 
cal switchesi 5 *^ or as gigahertz oscillatorsit 

In this paper, we study theoretically nanoelectrome- 
chanical effects in doubly-clamped suspended carbon 
nanotubes. Doubly-clamped suspended single- and 
multi-wall carbon nanotubes have been previously fab- 
ricated, and their transpor1pi2ii2i22i2i, acoustoelectrio22, 
thermal^, and elastio24 properties have been measured. 
We consider a single-wall carbon nanotube (SWNT) in 
which Coulomb-blockade effects dominate transport, and 
demonstrate that a gate manipulates the tube in an ef- 
fective way. The applied gate voltage bends the tube, 
changes the stress and thus influences the electric and 
mechanical properties. 

This paper is organized as follows: The next Section 
describes the model with inclusion of the influence of 
initial stress and thermal fluctuations. We concentrate 
on the case where the junction capacitances are zero so 
that analytical expressions are obtained. Section [ffi] de- 
scribes the influence of nanoelectromechanical effects on 
Coulomb blockade and shows that intrinsic bistability 
occurs when the tube is strained. Section IIVI discusses 
the eigenmodes and the influence on the initial strain on 
them. In Section Ivl junction capacitances are no longer 
neglected and we also show the effect of a non-uniform 
charge distribution. We end with some remarks on the 
limitations of our model. 



FIG. 1: A schematic drawing of a suspended nanotube capac- 
itively coupled to a gate and clamped on both sides to metal 
pads that serve as tunnel contacts. A voltage V is applied to 
the left pad. 

Studies with NEMS have mostly been performed in 



II. DISPLACEMENT, STRESS, AND ENERGY 

A. Equilibrium position 

We consider a SWNT (modeled as a rod of length L 
along the x-axis), freely suspended between source and 
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drain electrodes, in the vicinity of a gate (see Fig. 1). 
The nanotube is attached to the electrodes via tunneling 
contacts. An electrostatic force (gate voltage) bends the 
tube; the deviation from a straight line is denoted by z(x) 
with < x < L. The elastic energy of the bent tube is— 
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where E, I = irr /4, and 5* = irr 2 are the elastic mod- 
ulus, the inertia moment and the cross-section, respec- 
tively. Here, r is the (external) radius of the tube. The 
first term in Eq. is the energy of an unstressed bent 
rod; the two other terms describe the effect of the stress 
force T = Tq + T. Here To is the residual stress which 
may result e.g. from the fabrication, and the induced 
stress T is due to the elongation of the tube caused by 
the gate voltage, 
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To write down the electrostatic energy, we denote the 
capacitances of the barriers connecting the nanotube 
with the source and drain as Cl and Cr, respectively 
(see Fig. 1). The capacitance to the gate per unit length 
is c(z). Approximating the gate by an infinite plane at a 
distance R from the nanotube, we obtain 
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where the Taylor expansion restricts validity to z <C R. 
In this limit van der Waals forces between the nanotube 
and the substrate can be neglected. The electrostatic 
energy of the system reads 



W est [z(x)} = 
C L (C R 
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where V and V G are the potentials of the source and the 
gate respectively (the drain potential is set to zero) , ne is 
the (quantized) excess charge on the nanotube, and for a 
uniform charge distribution the capacitance to the gate 
equals 



C G = 



z(x)]da 



Note, that the last term in Eq. depends on the 

tube displacement and thus on the number of electrons. 
Therefore, it can not be omitted as in the standard 
Coulomb blockade treatment that replaces this term by a 
constant making W es t a periodic function of gate voltage. 



In the following, we concentrate on the analytically 
tractable case C^, Cr = 0. The general case is considered 
in Section For a moment, we also assume To = 0. In 
this situation, the expression for the electrostatic energy 
simplifies, 
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Minimizing the energy, 

W n [z(x)] = W e l[z(x)] + W est [z(x)], 

with respect to z, one finds the equation determining the 
tube position 2 ^, 
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where Kq is the electrostatic force per unit length, which 
we approximate by a constant. Higher-order terms are 
small for z <C R. To solve Eq. © we have to assume 
that the stress force T is constant, and find it later from 
the self-consistent condition (Eq. @). 

The solution of Eq. © with the appropriate boundary 
conditions (for the doubly-clamped rod z(0) = z(L) = 
z'(0) = z'(L) = 0) has the form 
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Substituting this into Eq. @, a relation between the 
stress T and the external force Kq is obtained. In the 
limiting cases, it reads 
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The first line corresponds to weak bending of the tube: 
The energy associated with the bending exceeds the en- 
ergy of the stress. Generally, it is realized for z < r. 
The second line describes strong bending, when the tube 
displacement is large (r < z R,L). 

For the displacement of the tube center z™ ax — 
z n (L/2) we find 
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For a SWNT with r = 0.65 nm, E = 1.25 TPa, L = 
500 nm and R = 100 nm (to be referred to as the E- 
nanotube) the crossover from weak to strong bending, 
T ~ EI/L 2 , occurs already at n ~ 5 -j- 10. In the strong- 
bending regime, the displacement of the E-nanotube is 
(in nanometers) z™ ax = 0.24n 2//3 . Note that this regime 
is not accessible with state-of the art silicon submicron 
devices, which are always in the weak-bending limit. 
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FIG. 2: Calculated displacement as a function of gate voltage 
for the E-nanotube: r = 0.65 nm, E = 1.25 TPa, L = 500 nm 
and R = 100 nm. At Vg ~ 0.5V , there is a crossover from 
weak bending with a Vj-dependence to strong bending with 
a Vq 3 dependence. 



B. Charge and energy 

For comparison with experiments, we have to relate the 
charge ne to the gate voltage by minimizing the energy. 
The expression for the energy (elastic plus electrostatic) 
of the tube at equilibrium in the limiting cases reads 
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T < EI/L 2 ; 
T> EI/L 2 . 



The first two terms represent the electrostatic energy of 
a straight tube, and the third one is due to the elastic 
degrees of freedom (stress, bending, and change of Cq 
due to displacement). This nonlinear, nanomechanical 
term is typically a small correction: For the E-nanotube 
it becomes of the same order as W st if n ~ 3000 in which 
case Eq. © is not valid anymore. The negative sign 
of the nanomechanical contribution is easily understood: 
As the gate voltage changes, the movable tube adjusts 
not only its charge, but also its position, which leads to 
a lower energy as compared to the fixed-position system. 
The value of n which minimizes the energy is 
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with Int denoting the integer part of the expression. The 
small correction Sn in the strong-bending regime is pro- 

5 /3 

portional to V G . Thus, the tube displacement z max 
changes in discrete steps when Vq is varied as shown in 
Fig. 2. The envelope is proportional to Vj (weak bend- 



ing) or Vq 3 (strong bending). In the absence of charg- 
ing effects and tension, the displacement is given by the 
dashed line as previously found in simulations of Ref. ITU 



The preceding considerations are restricted to the case 
of zero temperature. To understand the role of the tem- 
perature, we now evaluate the effect of thermal fluctua- 
tions on the equilibrium position of the tube. 

The variance of the position of the tube center at a 
given charge n can be generally represented as a func- 
tional integral, 

varz„ = ([z(L/2)-z n (L/2)} 2 ) 

= gj2 D *( X ) eX P [-W n [z]/k B @ + Jz{L/2)\ 

1 -1 

Dz(x)exp(-W n [z}/k B Q) , (11) 

where & is the temperature. Except for n = 0, the func- 
tional integral in Eq. I|ll|) is not Gaussian and has to 
be linearized around the equilibrium solution z n (x), Eq. 
0. The remaining Gaussian integral can be calculated, 
and we arrive at 

var z n = fc B 9C(i/2), (12) 
where £(x) solves the equation 

EK"" - |f / z' 2 dx C" - ^-4 J ('z> n dx 

= S(x-L/2). (13) 

In the two limiting cases of weak and strong bending, 
the solution of Eq. (|13f) yields 
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where the stress T is still given by the lower line of Eq. 
©. Thus, the fluctuations in the tube position are ex- 
pected to grow linearly with temperature. However, their 
magnitude is small. For the E-nanotube, at 100K the 
fluctuations in the n = state are of the order of 0.1 
nanometer, and at least an order of magnitude less in 
the strong-bending regime. 

In the calculations, we have assumed that the charge 
ne is a fixed quantity. Close to the degeneracy points 
W£ q = W^+i thermal fluctuations may induce switching 
between the states with charges ne and (n+ l)e, in which 
case Eq. (11411 is no longer valid. However, the range of 
voltages where switching is important, is narrow. 



III. COULOMB EFFECTS AND BISTABILITY 

A. Coulomb blockade 

Since the nanotube is attached to the electrodes by 
tunneling contacts, it is in the Coulomb blockade regime. 
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We define the energy to add the nth electron to the tube 
as S n = W n — W n -\. Then, if the nanotube contains 
n > electrons, the conditions that current can not flow 
(is Coulomb blocked) are S n < 0,eV < S n +i- In quan- 
tum dots, S n depends linearly on the bias V and gate V g 
voltages. Thus, in the Vg — V plane regions with zero 
current are confined within Coulomb diamonds, that are 
identical diamond-shape structures repeating along the 
Vg-bxis. 

In a suspended carbon nanotube, in addition to the 
purely Coulomb energy, we also have the nanomechanical 
corrections. Generally, these corrections make the rela- 
tions between V and Vg, which describe the boundaries 
of Coulomb blockade regions, non-linear. Consequently, 
the Coulomb "diamonds" in suspended nanotubes are 
not diamonds any more, but instead have a curvilin- 
ear shape (with the exception of the case Cl = Cr = 
0). Their size is also not the same and decreases with 
\Vg\- Thus, the mechanical degrees of freedom affect the 
Coulomb blockade diamonds. However, since these ef- 
fects originate from the nanomechanical term which is 
typically a small correction, its influence on Coulomb di- 
amonds is small as well. For the E-nanotube, these effects 
do not exceed several percents for typical gate voltages. 



B. Two-gate setup and bistability 

To demonstrate that the nanomechanical effects can 
not generally be omitted, we consider a suspended tube 
symmetrically placed in between two gates and show be- 
low that bistability in the tube position occurs^. 

Fig. 1 again presents the schematic setup, but the sus- 
pended tube is placed between two gates, labeled up (U) 
and down (D). Since up and down capacitances are con- 
nected in parallel, their sum Co — Cjj + Cd matters. 
Assuming that the distance of the straight tube to both 
gates is the same, we write 
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Expanding this for z <C R and calculating the electro- 
static force, we arrive at an equation similar to Eq. JfjJ), 
with a constant force Kq that is replaced by jz, where 
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We now solve this equation in the strong-bending regime. 
For this purpose-- we disregard the term IEz"", and 
use the boundary conditions z(Q) = z(L) — 0. Multiple 
solutions emerge; the ones with the lowest energy are 
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Thus, the tube in the strong-bending regime can oscil- 
late between the two symmetric positions. This creates 



a basis for observation of quantum effects, as discussed 
in Ref. |2g. We emphasize once again, that within this 
model, the multi-stability is due to the charging of the 
tube in combination with the non-linearity. 



IV. EIGENMODES 

The eigenfrequency of a particular eigenmode is an 
important directly measurable^ property. In future ex- 
periments on suspended tubes we expect that the eigen- 
modes influence tunneling ("phonon-assisted tunneling") 
in a similar way as observed for a single Cqq molecule^. 
Below, we demonstrate that the effect of the electrostatic 
interactions on the elastic properties (specifically, eigen- 
frequencies) is strong and changes the behavior qualita- 
tively. 

To find the eigenmodes, we apply a gate voltage with 
a large dc (single gate) and a small ac component. The 
displacement z(x, t) is time-dependent, which provides 
an external force —pSz to Eq. ©, where p equals 1.35 
g/cm 3 . Eq. ijfJJl must be solved first with a constant 
stress, and then the stress is found self-consistently. The 
tube displacement has a small ac component Sz on top 
of a large static one. The self-consistency procedure is 
essentially the same and again leads to Eq. JHJ. Thus, the 
dc component of the gate voltage determines the stress 
T and it therefore controls the eigenmodes. 

The frequencies of the (transverse) eigenmodes are 
found from the requirement that the equation 



pSco 2 Sz = 



(17) 



with the boundary condition Sz(0) = Sz(L) — Sz'(0) = 
Sz'(L) = has a non-zero solution. This yields the fol- 
lowing equation for the frequency to, 



1 y 2 — yi, 
cosh yi cos y 2 - 
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In the following, we restrict ourselves to the fundamental 
(lowest frequency) eigenmode u>o. In the limiting cases, 
the solutions of Eq. (|18(l are 



(19) 



The second terms on the rhs represent small corrections 
to the first ones. 

The frequency dependence luo oc L~ 2 is associated with 
a loose string, while ujq cx L _1 means that the string is 
tied like in a guitar. Our results show that the behavior 
of the tube crosses over from "loose" to "tied" as Vg in- 
creases. For the fundamental mode, the crossover occurs 
at £L ~ 1, corresponding to the crossover from weak 
to strong bending. The middle curve in Fig. 3 shows 
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the frequency of the fundamental mode as a function of 
gate voltage (zero residual stress). The arrow denotes 
the cross-over from weak to strong bending. 

The gate voltage dependence of the frequency is a step- 
wise function, as shown in the inset of Fig. 3. Steps 
occur whenever an additional electron tunnels onto the 
tube. For the E-nanotube, their height is ~5 MHz, which 
is measurable. Note, that the present submicron silicon 
devices are always in the weak-bending regime so that 
corrections due to the second term in Eq. (|19J) are too 
small to be measured. Furthermore, one should realize 
that frequency quantization is only observable if the fre- 
quency itself is greater than the inverse tunneling time 
for electrons. 

We now consider the effect of a residual stress (To ^ 0). 
First, we obtain the stress by solving Eqs. (2), © (in 
the latter, T is replaced by T + To). In particular, for 

a negative stress T + T < 0, T EI/L 2 , Eq. © 

acquires several solutions. This signals Euler instability: 
the tube bends in the absence of an external force. 



It further drops to zero at the Euler instability threshold. 

The qualitative difference between the various regimes 
means that by measuring the gate voltage dependence of 
too one can determine the sign of To and get a quantitative 
estimate. On the other side, the gate effect can be used 
to tune the eigenfrequencies. We also mention that in 
the absence of charging effects, the steps vanish but the 
overall shape of the curves in Fig. 3 remains the same. 

V. RELAXING THE APPROXIMATIONS 

While considering equilibrium displacement and eigen- 
modes of the nanotube, we made a number of simplifying 
approximations. In this Section, we consider two of them 
— disregarding the capacitances Cl,r and uniform dis- 
tribution of the charge — and show that relaxing these 
approximations affects the above results quantitatively, 
but not qualitatively. 

In this Section, we consider the case of zero residual 
stress T = 0. 



A. Finite capacitances to the leads 
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FIG. 3: Gate voltage dependence of the frequency ujq of the 
fundamental mode for three different values of the residual 
stress. Numbers are taken for the E-nanotube (see Fig. 2). 
The fundamental mode of an unstressed tube is 140 MHz (thin 
horizontal line). The inset is an enlargement of the To = 
curve of the main figure showing step-wise increases of ujq 
whenever an additional electron tunnels onto the tube. 



We now relax the limitation Cl, Cr = 0. For the gen- 
eral case, Eq. (|HJ| still holds, however, the force Kq must 
be adjusted, 
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(20) 



where Co = L/(2\n2R/r) is the capacitance of the 
straight nanotube to the gate. The results of the nu- 
merical solutions for the displacement and the frequency 
of the fundamental mode are plotted in Fig. 0] For sim- 
plicity we have taken Cl = Cr = §Cq\ the four curves 
correspond to different values of the parameter <f>. The 
curves with = are the same as the ones in Figs. |5J |3] 
The plots demonstrate that the qualitative picture re- 
mains the same if we include finite capacitances to the 
leads. The steps observed for = become skewed with 
the increase of Cl and Cr (see inset of Fig. 4). At a 
certain they disappear. For > 10 the plots are, on 
the scale presented, the same. 



If the residual stress is large, To ^ EI/L 2 , the tube al- 
ways acts like a tied string (upper curve in Fig. 3). The 
frequency depends weakly on Vg for low voltages, and 
above T ~ To (denoted with the arrow) grows with an 

2/3 

envelope oc V G . For negative To the picture is qualita- 
tively different (lower curve in Fig. 3). Whereas for large 

2/3 

gate voltages the envelope is still proportional to V G , 
the frequency dives below the value for an unstressed 
tube (22.38(EI/pS) 1/2 L~ 2 , represented by the thin solid 
line in Fig. 3), when the overall stress becomes negative. 



B. Non-uniform charge distribution 

Above, we have assumed a uniform charge distribution 
along the nanotube. Rather than trying to analyze the 
effect in general, we consider the opposite situation when 
the excess charge is concentrated at one point (to be more 
precise, in a concise region of the tube radius r), which 
may represent, for instance, a pinning center. This center 
is placed in the middle of the nanotube. Though we be- 
lieve that the charge distribution in suspended nanotubes 
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FIG. 4: Above: Displacement as a function of gate voltage 
for the E-nanotube with finite capacitances to the leads. The 
four curves correspond to different values of the parameter <f>, 
defined as Cl = Cr = 4>Cg- The inset is an enlargement 
of the main figure. Below: The frequency of the fundamen- 
tal mode normalized to the fundamental frequency of an un- 
stressed tube n = (22.38L- 2 (EI/pS) 1/2 = 141 MHz for the 
same parameters as above. 



is closer to uniform, this situation applies to a suspended 
quantum dot as realized recentlySi. 

The gate-charge capacitance Cq in this geometry is 
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1 



i i_ ' 

r 2R 



and we proceed to obtain the equations of motion, 

IEz»»-Tz» = F5(x-±\, F = 

\ 2 J ' 4R 2 ' 



(21) 



(22) 



where we again set Cl = Cr = 0. 

The solution with the same boundary conditions as 
previously, z(0) = z(L) = z'(Q) = z'{L) = 0, and with z, 
z' , and z" all continuous at x = L/2, is 



z(x) 



F 



(23) 



\ <j tanh — [cosh t;x — 1] — sinh £x + t;x 



for < x < L/2. For L/2 < x < L the coordinate x 
should be replaced by (L — x) because z(x) — z(L — x). 
As before, £ = {T/EI) 1 / 2 and Eq. @ is used to obtain 



FIG. 5: Displacement of the center of the E-nanotube (above) 
and the frequency of the fundamental eigenmode (below) for 
uniform charge distribution and for the case that the charge 
is concentrated at one point. Q. is the fundamental frequency 
of an unstressed tube. 



the stress self-consistently, 



T = 



F 2 L 4 S/(30720EI 2 ), T < EI/L 2 , 
(1/2)(£SF 2 ) 1/3 , T^EI/L 2 . 



(24) 



Consider now the strong-bending regime and compare 
the results for the stress T u for the uniform (lower line 
of Eq. JHJ and T n for the concentrated (lower line of Eq. 
(|24|l charge distributions, 



T — T 



V3L 
4R 



2/3 



(25) 



For L 3> R we formally have T n 3> T u . This means 
that for the same gate voltage more stress is induced at 
the nanotube if the charge is concentrated at one point. 
Also, the displacement of the tube is greater in the con- 
centrated case, 



0.87 



1/3 



Thus, if the charge distribution is concentrated, NEMS 
are "more effective" than for the uniform charge. For 
the E-nanotube the ratio between non-uniform and uni- 
form maximal displacement is 1.49. The difference be- 
tween uniform and non-uniform charge distributions is 
illustrated in Fig. [S] 
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VI. DISCUSSION 

The presented model is simplified in many respects. 
Mechanical degrees of freedom are introduced via clas- 
sical theory of elasticity: The nanotube (modeled by a 
rod) is considered as incompressible and without internal 
structure. This is justified, since so far the theory of elas- 
ticity has described all existing experiments on carbon 
nanotubes well. For SWNTs it has been also supported 
by simulations (Ref. \13t) . Creation of defects in SWNT 
starts at deformations of order of ten percents. For larger 
deformations (see e.g. Ref. Esh we expect strong devi- 
ations from the behavior we describe, but this typically 
lies outside our applicability range z <C R- We have 
neglected damping, which is also expected to originate 
from the creation of the defects and to be irrelevant in 
this range. We also disregarded quantum effects (cotun- 
neling and finite spacing of quantum levels of electrons 
in the tube). These issues need to be clarified for a de- 
tailed comparison with the experimental data, and will 
be a subject for future research. 

Our main result is that the nanotube can be manipu- 



lated by the gate voltage, which determines its deforma- 
tion and stress, and modifies the eigenmodes. Though 
the eigenmodes of nanotube ropes have been measured 
in Ref. three years ago, the strain dependence of the 
eigenmodes was only recently reported in Ref. |29| which 
was published after this manuscript had been submit- 
ted for publication. Ref. |22l demonstrates this effect for 
singly-clamped multi-wall carbon nanotubes. We expect 
that our predictions will soon be tested in experiments 
on doubly-clamped SWNTs. 

We mention also one more paper published after the 
submission of our manuscript, Ref. 1271 which shows mea- 
surements on a suspended quantum dot. Though the 
focus of our study was on carbon nanotubes, all the cal- 
culations can be immediately applied to this case as well. 

We thank Yu. V. Nazarov, P. Jarillo-Herrero, 
L. P. Kouwcnhoven and C. Dekker for discussions. This 
work was supported by the Netherlands Foundation for 
Fundamental Research on Matter (FOM) and ERATO. 
HSJvdZ was supported by the Dutch Royal Academy of 
Arts and Sciences (KNAW). 
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